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ABSTRACT 

How massive stars die - what sort of explosion and remnant each produces - depends chiefly 
on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar 
winds are the only means of mass loss, and these are chiefly a function of the metallicity of the 
star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects 
the evolution and final fate of massive stars. We map, as a function of mass and metallicity, 
where black holes and neutron stars are likely to form and where different types of supernovae 
are produced. Integrating over an initial mass function, we derive the relative populations as a 
function of metallicity. Provided single stars rotate rapidly enough at death, we speculate upon 
stellar populations that might produce gamma-ray bursts and jet-driven supernovae. 

Subject headings: massive stars, supernovae, stellar remnants, neutron stars, black holes, gamma-ray 
bursts, collapsars 

1. Introduction history of its mass loss. For single stars, mass loss 

occurs as a result of stellar winds for which there 
The fate of a massive star is governed chiefly g^ist semi-empirical estimates. Thus, within cur- 

by its mass and composition at birth and by the ^g^tiy existing paradigms for the explosion, the 
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fate of a star of given initial mass and compo- 
sition is determined (the Russell- Vogt theorem). 
If so, one can calculate realization frequencies for 
stellar explosions and remnants of various kinds 
and estimate how these might have evolved with 
time. 

Such estimates are fraught with uncertainty. 
The litany of complications is long and requires 
discussion (§ 6). No two groups presently agree, 
in detail, on the final evolution of any massive star 
(including its explosion energy, remnant mass, and 
rotation rate) and the scaling of mass loss with 
metallicity during different evolutionary stages is 
widely debated. Still it is worthwhile to attempt 
an approximate table of histories. We would like 
to know, within the comparatively well under- 
stood domain of stars that do not experience mass 
exchange with a companion and for a particu- 
lar set of assumptions regarding mass loss and 
explosion, what sort of supernova each star pro- 
duces and what sort of bound remnant, if any, it 
leaves. If possible, we would also like some indi- 
cation of which massive stars might make gamma- 
ray bursts. 

In this paper we construct such a table of stel- 
lar fates and remnants. In § 2 we describe our 
assumptions regarding mass loss, explosion mech- 
anism(s), and remnant properties, and in § 6 dis- 
cuss the uncertainties. Section 4 delineates the 
sorts of stellar explosions and collapses we want 
to distinguish, and in § 5 we discuss the resulting 
realizations of different outcomes as a function of 
metallicity in the galaxy. 

2. Assumptions 

2.1. Stellar Models and Paradigms 

The stellar models used in this paper were taken 
from Woosley & Weaver (1995); Heger & Woosley 
(2002b); Woosley et al. (2002) and Heger et al. 
(2003a). These papers treat the evolution of mas- 
sive stars in the range 9 to 300 M© calculated with- 
out rotation from birth on the main sequence to 
death, either as iron-core collapse supernovae (he- 
lium core masses at death less than about 65 Mq) 
or pair instability supernovae (helium core masses 
at death greater than 65 Mq and up to about 135 
Mq). The effects of mass loss were included in 
those studies as discussed in § 2.2. 

We shall presume here that the explosion mech- 



anism, however it may operate, and the remnant 
properties are determined by the mass of the he- 
lium core when the star dies. (Perhaps the carbon 
oxygen core mass is a better discriminant, but sys- 
tematics of the two are very similar). As the mass 
of the hehum core increases, so does its binding 
energy and entropy. Because of its higher entropy, 
a larger helium core also has, on the average, a 
larger iron core mass, and a shallower density gra- 
dient around that core (Woosley et al. 2002). Con- 
sequently such stars are harder to explode (Fryer 
1999; Fryer & Kalogera 2001). Even in "success- 
ful" explosions where a strong outward shock is 
born, mass may later fall back onto a neutron star 
remnant turning it, within one day, into a black 
hole. We thus distinguish black holes that are pro- 
duced promptly or "directly" from those made by 
fall back. 

Fryer (1999) has estimated that the helium core 
mass where black hole formation by fall back en- 
sues is about 8 M© (a <25 M© main sequence star) 
and that direct black hole formation occurs for he- 
lium cores over 15 M© (40 M© main sequence star 
with no mass loss). These numbers are uncertain 
(§ 6.2), but are representative choices. It is as- 
sumed that a baryonic remnant mass of over 2.0 
M© will produce a black hole. 

While the helium core mass governs the ex- 
plosion mechanism, the hydrogen envelope is 
largely responsible for determining the spectrum 
(at peak) and light curve of common Type II su- 
pernovae. Stars with massive hydrogen envelopes 
when they die will be Type Hp; low mass en- 
velopes will give Type IIL and lib; etc. (§ 4). 
An exception are supernovae of Types lb and Ic 
whose light curves do depend sensitively on the 
helium core mass since all the hydrogen envelope 
has been removed. The light curves of Types lib, 
lb, Ic, and 87A-like explosions are also sensitive 
to the amount of ^^Ni made in the explosion. 

2.2. Mass loss 

The principal physics connecting the final evo- 
lution of a star to its metallicity is its mass loss. 
Low metallicity stars have less mass loss and have 
bigger helium cores and hydrogen envelopes when 
they die. To a lesser extent, metallicity also affects 
whether the presupernova star is a red or blue su- 
pergiant (Langer & Maeder 1995). 
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For main sequence stars and red supergiants 
the mass loss rates employed in the studies cited 
above were taken from Nicuwenhuijzen & de Jager 
(1990). For Wolf-Rayet stars, a mass-dependent 
mass loss rate (Langer 1989) was assumed using 
the scaling law established by Brown (1997); Well- 
stein & Langer (1999), but lowered by a factor 3 
(Hamann & Koesterke 1998). Wind-driven mass 
loss is believed to be metallicity dependent and a 
scaling law oc \fZ has been suggested for hot stars 
(Kudritzki 2000; Nugis & Lamers 2000). Woosley 
et al. (2002) assumed that the same scaling law 
holds for Wolf-Rayet stars (Vanbeveren 2002) and 
blue and red supergiants as well. "Metallicity" is 
assumed here to be the initial abundance of heavy 
elements, especially of iron, not the abundances of 
new heavy elements like carbon and oxygen in the 
atmospheres of WC and WO stars (§ 6.1). 

In very massive stars above '--'60 Mq, the cp- 
silon mechanism for pulsational driven mass loss 
sets in and enhances the mass loss during central 
hydrogen burning. Opacity-driven pulsations also 
become important, if not dominant, at high metal- 
licity (Baraffe et al. 2001). At very low metallic- 
ity on the other hand, Baraffe et al. (2001) have 
shown that primordial stars should not have sig- 
nificant mass loss due to pulsations. This suggests 
significant evolution in the mass loss of very mas- 
sive stars with metallicity (Figs. 1-4). 

3. Remnant Properties 

Fig. 1 shows the expected remnant types as a 
function of mass and initial "metallicity" for the 
above assumptions. In preparing Fig. 1, it is as- 
sumed that stars below ~ 9 M0 do not form mas- 
sive enough cores to collapse, that they end their 
lives as white dwarfs. Just above this mass lies 
a narrow range, 9 — 10 Mq, where degener- 
ate oxygen-neon cores are formed that cither col- 
lapse due to electron capture (Barkat et al. 1974; 
Miyaji et al. 1980; Nomoto 1984; Habets 1986; 
Miyaji & Nomoto 1987; Nomoto 1987; Nomoto & 
Hashimoto 1988) and make a neutron star or lose 
their envelopes and make white dwarfs (Garcia- 
Berro & Iben 1994; Ritossa et al. 1996; Garcia- 
Berro et al. 1997; Iben et al. 1997; Ritossa et al. 
1999). Above ^ 10 M0 core collapse is the only 
alternative. 

Wherever this transition between white dwarf 



formation and iron core collapse lies, it should de- 
pend very little on metallicity and thus appears 
as a vertical line in Fig. 1. At low metallicities, 
the boundaries for black hole formation are also 
defined entirely by the initial stellar mass since 
there is a one to one correspondence between ini- 
tial stellar mass and final helium core mass. 

For stars of higher metallicity, mass loss be- 
comes increasingly important resulting in smaller 
helium cores for a given initial mass. If the 
star loses its entire hydrogen envelope (to the 
right of the green line in Figs. 1-4), its rate of 
mass loss increases significantly (e.g., Langer 1989; 
Hamann et al. 1995) producing much smaller he- 
lium cores at collapse. This effect underlies the 
abrupt change in the otherwise vertical bound- 
aries between neutron star, fallback black hole 
and direct black hole formation. For very massive 
stars, the remnant of the collapsing star depends 
sensitively on the metallicity. Above 40 M0, low 
metallicity stars form black holes directly, while 
at higher metallicities black holes of smaller mass 
are produced by fall back until, ultimately, only 
neutron stars are made. Winds are assumed to be 
stronger in higher mass stars, so the metallicity at 
which these transitions occur decreases with mass. 
But beyond ^ 100 Mq, this limit may rise again 
due to high enough initial mass or the significant 
role of evolution phases with lower mass loss rates 
(e.g., a WNL phases; see Brown et al. 2001). 

At low metallicities, there is also a range of 
masses for massive stars that leave behind no rem- 
nant whatsoever. These are the pair-instability 
supcrnovae. If the helium core exceeds ~ 65 M0, 
corresponding to a ~ 140 Mq initial mass for stars 
without mass loss, the pulsational pair instabil- 
ity (Heger & Woosley 2002b) becomes so violent 
that the star is disrupted entirely. When the he- 
lium core mass at the end of central carbon burn- 
ing exceeds ~ 135 Mq for non-rotating stars (ini- 
tial mass of ~ 260 Mq without mass loss), photo- 
disintegration in the center leads to collapse to a 
very massive black hole (> 100 Mq), once again 
forming a black hole directly (Fryer et al. 2001; 
Heger & Woosley 2002b). However, as the metal- 
licity increases, mass loss shifts the regime of pair- 
instability supernovae to higher initial masses. At 
still higher metallicities, these supcrnovae do not 
occur at all (Baraffe et al. 2001) because the pro- 
genitor stars are pulsationally unstable. 
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4. Supernovae 

4.1. Supernovae of Type Hp and IIL 

It has long been reeognized that massive stars 
produce supernovae (Baade & Zwicky 1938). In 
this paper, we assume the foUowing progenitor 
properties for the different core collapse supernova 
types: 



SN Type 


pre-SN stellar structure 


Hp 


> 2Mo H envelope 


IIL 


< 2M0 H envelope 


Ib/c... 


no H envelope 



The lower and upper limits of main sequence 
mass that will produce a successful supernova 
( "M-lower" and "M- upper" ) — one with a strong 
outgoing shock still intact at the surface of the star 
— has long been debated. On the lower end, the 
limit is set by the heaviest star that will eject its 
envelope quiescently and produce a white dwarf. 
Estimates range from 6 to 1 1 M0 with smaller 
values characteristic of calculations that employ 
with a large amount of convective overshoot mix- 
ing (Marigo et al. 1996; Chiosi 2000) and the upper 
limit determined by whether helium shell flashes 
can eject the envelope surrounding a neon-oxygen 
core in the same way they do for carbon-oxygen 
cores (§ 3). It may also slightly depend on metal- 
licity (Cassisi & Castellani 1993). Here we will 
adopt 9 Mq for M-lower. 

The value of M-upper depends on details of the 
explosion mechanism and is even more uncertain 
(§ 6.2). Fryer & Kalogera (2001) estimate 40 Mq, 
but calculations of explosion even in supernovae 
as light as 15 M0 give widely varying results. It is 
likely that stars up to at least 25 Mq do explode, 
by one means or another, in order that the heavy 
elements be produced in solar proportions. The 
number of stars between 25 and 40 Mq is not large. 
Here we have taken what some may regard as a 
rather large value. M-upper equals 40 Mq (Fig. 2). 

For increasing metallicity mass loss reduces the 
hydrogen envelope at the time of core collapse. 
A small hydrogen envelope (< 2Mq) can't sus- 
tain a long plateau phase in the light curve, and 



only Type IIL supernovae or, for very thin hydro- 
gen layers. Type lib supernovae result (Barbon 
et al. 1979; Filippenko 1997). It is also necessary 
for Type IIL supernovae that the radius be large 
(Swartz et al. 1991) and helpful if the ^^Ni mass is 
not too small. The minimum metallicity for Type 
IIL supernovae in single stars, is set by the re- 
quirement that the mass loss needs to be strong 
enough to remove enough of the hydrogen enve- 
lope (Fig. 2). In single stars Type IIL/b SNe are 
formed only in a thin strip where the hydrogen 
envelope is almost but not entirely lost. Gaskell 
(1992) finds that Type IIL supernovae are cur- 
rently about 10% - 20% as frequent as Type Hp. 

For increasing metallicity this domain shifts 
to lower initial mass. Below a certain minimum 
metallicity we do not expect Type IIL supernovae 
from single stars at all. Indeed, those stars that 
form at the lowest (possible) metallicities will be 
so massive that they frequently form black holes 
by fall back and have not very luminous super- 
novae. This will be particularly true if the stars 
explode as blue supergiants but lack radioactivity. 

4.2. Type lb and Ic Supernovae 

A complication is that Type Ib/c SNe with 
masses above 4-5 Mq , which may be the most com- 
mon ones to come from single stars, also have dim 
displays even if they are still powerful explosions 
(Ensman & Woosley 1988), i.e., the progenitor 
stars' cores are not so massive that they encounter 
significant fallback. In this paper, we do not dif- 
ferentiate these types of supernovae from our set 
of normal supernovae. Our assumptions regarding 
the different types of supernovae are summarized 
in the table below: 



Type Ib/c 
He core mass 
at explosion 


explosion 
energy 


display 


> 15 Mq 


direct collapse 


none^ 


- 15-8Mq 


weak 


dimt 


- 8 - 5 Mq 


strong 


possibly dim 


< 5Mq 


strong 


bright 



^if not rotating 
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Clearly, mass loss is a key parameter and both 
high metallicities (and high initial masses) are re- 
quired to produce Type Ib/c supcrnovae in sin- 
gle stars. Woosley et al. (2002) find that for so- 
lar metallicity the limit for non-rotating stars is 
~ 34 M0 . These supernovae can be weak and their 
later fallback will produce BH remnants. As with 
the Type II black-hole forming supernovae, we an- 
ticipate that this fallback, in particular of the ^^Ni 
lost this way, may weaken the brightness of the su- 
pernova display, similar to the case of weak Type 
II SNc. 

4.3. Nickel-deficient Supernovae 

The light curve of most supernovae is a con- 
sequence of two energy sources - shock-deposited 
energy and radioactivity, especially the decay of 
5^Ni to ^^Fe. There are cases however, where 
the radioactive component may be weak or ab- 
sent. If the hydrogen envelope is still present, a 
bright supernova may still result with the bright- 
ness depending on the explosion energy (Popov 
1993), but the light curve lacks the characteris- 
tic radioactive "tail" (e.g., SoUerman et al. 1998; 
Turatto et al. 1998; Benetti et al. 2002). If the hy- 
drogen envelope is is gone (Type Ib/c), the con- 
sequences for the light curve are more dramatic 
and the supernova may be for practical purposes 
invisible. 

Four cases of nickel-deficient supernovae may 
be noted. 

1) Stars in the mass range 9 to 11 Mq. Such 
stars have steep density gradients at the 
edge of degenerate cores. The shock wave 
from core collapse heats very little material 
to greater than 5 x 10^ K and very little 
(<0.01Mq) ^^Ni is ejected (Mayle & Wilson 
1988) 

2) Stars which make ^^Ni but where the ^^Ni 
falls back into the remnant. This occurs 
for more massive stars with the threshold 
mass dependent upon both the presuper- 
nova structure and the explosion mechanism 
and energy. The boundary here is somewhat 
fuzzy because of the operation of mixing in 
conjunction with fallback. The lower limit 
for this regime is probably slightly larger 
than that for BH formation by fallback, the 



upper limit is where BHs are formed directly 
without initiating a supernova, i.e., 10 Mq < 
helium core mass < 15 Mq (stellar masses 
30 Mq < M < 40 Mq without mass loss). 

3) Pair-instability supernovae with helium core 
masses in the range 65 to < 85 Mq. Pair- 
instability supernovae, which probably only 
existed in the early universe can have light 
curves ranging from very faint if they have 
lost their hydrogen envelopes and eject no 
^^Ni to exceptionally brilliant if the converse 
is true (helium core > 100 Mq; Heger & 
Woosley 2002b; Heger et al. 2002a). 

4) Pulsational pair-instability supernovae with 
helium core masses in the range to 
65 Mq . This instability occurs after cen- 
tral carbon burning but before the collapse. 
Though each pulse can have up to several 
10^^ erg, only the outer layers of the star are 
expelled and contain no '^^Ni (see below). 

4.4. Pair-instability supernovae 

Very massive stars [M > 100 Mq) still form in 
the present galaxy (Najarro & Figer 1998; Eikcn- 
berry et al. 2001), but above ~ 60 Mq, nuclear- 
powered and opacity driven pulsations occur that 
increase the mass loss (e- and k- mechanisms) . Re- 
cently, Baraffe et al. (2001) have shown that both 
mechanisms are suppressed in extreme Pop III 
stars. Therefore it seems reasonable to assume 
that at sufficiently low metallicity [Z < 10~^Zq) 
very massive stars may retain most of their mass 
through the end of central helium burning, form- 
ing a massive helium core (Baraffe et al. 2001; Ku- 
dritzki 2002; Marigo et al. 2003). 

For zero-metallicity stars above ^ 100 Mq (he- 
lium cores > 42 Mq; Woosley 1986; Chiosi 2000; 
Heger & Woosley 2002b) stars encounter the the 
pair instability after central carbon burning (e.g.. 
Bond et al. 1984; Heger & Woosley 2002b). Be- 
tween 100 Mq and ^ 140 Mq (helium core mass 
< 65 Mq) the instability results in violent pulsa- 
tions but not complete disruption. The implo- 
sive burning is not energetic enough to explode 
the star. Depending on the mass of the star and 
the strength of the initial pulse, subsequent pulses 
follow after < 1 yr to > 10, 000 yr. These pulsa- 
tions continue until the star has lost so much mass, 
or decreased in central entropy, that it no longer 
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encounters the pair instability before forming an 
iron core in hydrostatic equilibrium. Since the iron 
core mass is large and the entropy high, such star 
probably finally make black holes. 

The typical energy of these pulses can reach a 
few 10^^ erg and easily expels the hydrogen en- 
velope, which is only loosely bound, in the first 
pulse (Heger & Woosley 2002b) - when these stars 
finally collapse they are thus hydrogen-free. Sub- 
sequent pulses may eject the outer layers of the 
helium core as well. Though the kinetic energy 
of these pulses may be well in excess of normal 
supernovae, they are less bright since they lack 
any ^^Ni or other radioactivities that could power 
an extended light curve. However, the collision of 
shells ejected by multiple pulses could lead to a 
bright display. 

For stars between ~ 140 and ~ 260 M0 (helium 
cores of ^ 64 to 133 Mq) the pair-instability 
is violent enough to completely disrupt the star 
in the first pulse (Ober et al. 1983; Bond et al. 
1984; Heger & Woosley 2002b). Explosion energies 
range from ~ 3x10^^ erg to < 10^^ erg and the 
ejected ^^Ni mass ranges from zero to > 50 Mq 
at the high-mass end (Heger & Woosley 2002b). 
Above ~ 260 Mq, the stars directly collapse to 
a black hole (Fryer et al. 2001; Heger & Woosley 
2002b) . Rotation would of course affect these mass 
limits. 

4.5. Very energetic and asymmetric su- 
pernovae 

4-. 5.1. Jet-powered supernovae 

A jet-driven supernova (JetSN) is a grossly 
asymmetric supernova in which most of the en- 
ergy comes from bipolar outflow from a central 
object. Though such supernovae may occur in as- 
sociation with gamma-ray bursts (GRBs), not all 
jet-powered supernovae will have sufficiently rela- 
tivistic ejecta to make such a hard display. The 
class of jet-powered supernovae is thus a broad one 
having GRB progenitors as a subset. 

Jet-Driven supernovae can be formed with or 
without hydrogen envelopes (MacFadyen et al. 
2001; Nomoto et al. 2002; Fig. 4). The hydrogen- 
free JetSNe are closely related to GRBs. Whether 
such stars produce JetSNe or GRBs (or both) de- 
pends upon the rotation and the explosion mech- 
anism. Until we understand both better, we can 



not distinguish between the two. 

4.5.2. Gamma-ray bursts and coUapsars 

The currently favored model for the forma- 
tion of gamma-ray bursts assumes that a narrowly 
beamed (6* < 10°) highly relativistic jet (F > 100) 
leaves a compact "engine" and produces 7-rays ei- 
ther by internal shocks or by running into some 
external medium (Frail et al. 2001). Currently two 
classes of GRBs are distinguished: long and short 
bursts (Fishman & Meegan 1995). It is assumed 
that the short class might originate from binary 
neutron stars (Eichler et al. 1989), the long class 
could be produced by the collapse of the core of 
a massive star (e.g., Popham et al. 1999). In the 
present work we adopt this assumption, focus on 
the long class of GRBs, and use "GRB" synony- 
mous for this class. 

The term "collapsar" is used to describe all 
massive stars whose cores collapse to black holes 
and which have sufficient angular momentum to 
form a disk. There are three possible varieties. 

I coUapsars that form black holes "directly" 
during the collapse of a massive core. Al- 
though the star collapses and initially forms 
a proto-neutron star, it is unable to launch a 
supernova shock and eventually (after ~1 s) 
collapses to form a black hole (Woosley 1993; 
MacFadyen & Woosley 1999). 

II coUapsars that form black holes by fallback 
after an initial supernova shock has been 
launched (MacFadyen et al. 2001). The ex- 
plosion is too weak to eject much of the 
star, and the subsequent fallback of mate- 
rial causes the neutron star in the core to 
collapse and form a black hole. 

HI coUapsars which do not form proto-neutron 
stars at all, but instead quickly collapse into 
massive black holes which grow through ac- 
cretion (Fryer et al. 2001). These coUapsars 
lead to the formation of massive ('^3OOM0) 
black holes. 



type 


time-scale 


energy 
budget 


initial 
BH mass 


I 


short 


low 


small 


II 


long 


low 


small 


HI 


long 


high 


large 
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The results can be summarized as: 

• Type I and II coUapsars without a hydrogen 
envelope can make ordinary GRBs, though 
those of Type II will tend to be longer. 

• Type II and III coUapsars without a hydro- 
gen envelope - maybe even with - can make 
very long GRBs (in their rest frame). 

• All three types can make bright jet-powered 
supernovae if a hydrogen envelope is present. 

Though we have described them as GRB pro- 
genitors, coUapsars probably produce a variety of 
outbursts from X-ray flashes to jet-driven super- 
novae. Calculations to reliably show which stars 
make GRBs as opposed to just black holes are 
presently lacking (though see Hcger & Woosley 
2002a). Here we will assume that coUapsars are 
made by some subset of those stars that make 
black holes (Fig. 3). 

It is agreed however, that coUapsars can only 
form GRBs if the star has lost its hydrogen enve- 
lope prior to collapse. Mass loss depends both on 
the stellar mass and metallicity and as both in- 
crease, the star uncovers more and more of its hy- 
drogen envelope. The green curve in Figs. 3 and 4 
denotes the boundary between stars which retain 
some of their hydrogen envelope and those that 
lose all of their hydrogen through mass loss. Above 

30 Mq, mass loss from winds become important, 
and as the initial mass of the star increases, lower 
and lower metallicities are required to retain the 
hydrogen envelope. Between 100 — 140 Mq, pulsa- 
tional instabilities are able to drive off the hydro- 
gen layers of the star, even at zero metallicities. 
This boundary, which determines where stars lose 
their hydrogen envelopes marks the lower bound 
for GRB producing coUapsars. The upper bound 
is set by those stars that collapse to form black 
holes. 

5. Stellar Populations 

With our evaluation of the possible fates of mas- 
sive stars from § 4, we estimate the distribution 
of compact remnants and of observable outbursts 
produced by these single stars. The results will 
be uncertain. Not only do the predictions depend 
sensitively on the regions outlined in Figs. 1-4, 



but also on the initial mas function (IMF) and its 
evolution. 

In Fig. 6, wc plot the fraction of massive stars 
forming neutron stars (solid line) and black holes 
[dotted line) assuming a Salpetcr IMF (Salpeter 
1955). At low metallicities, roughly 20% of mas- 
sive stars form black holes, and roughly 75 % 
form neutron stars. Half of those black holes 
form through fallback, the other half through 
direct collapse. Only 4 % of black holes form 
massive (>2OOM0) black holes. 1 % of massive 
stars form pair-instability supernovae (leaving be- 
hind no remnant whatsoever). As the metallic- 
ity increases, the fraction of stars producing black 
holes first increases slightly (as the pair-instability 
mechanism is shut off) and then decreases near so- 
lar metallicity as most massive stars lose so much 
mass that they collapse to form neutron stars in- 
stead of black holes. At these high metallicities, all 
black holes are formed through fallback. Note that 
direct collapse black holes are larger than fallback 
black holes and black holes will be larger, on av- 
erage, at low metallicity. In addition, if the black 
hole kick mechanism is powered by the supernova 
explosion, direct black holes will not receive kicks 
and these large black holes will tend to have small 
spatial velocities. 

There is increasing evidence that the IMF is 
more skewed toward massive stars (relative to a 
Salpetcr IMF) at low metallicities (e.g., Bromm 
et al. 2001; Abel et al. 2002, 2000). To include 
these effects, we have used the IMF for Popula- 
tion III by Nakamura & Umemura (2001). The 
thin lines in Fig. 6 show the change in the distri- 
bution of black holes and neutron stars using the 
Nakamura & Umemura (2001) IMF with the fol- 
lowing parameters: m^^ = 1.5, nip^ = 50, k = 0.5, 
a = P = 1.35 (see Nakamura & Umemura 2001 for 
details). We employ this IMF up to a metallicity 
that corresponds to the last occurrence of (non- 
pulsational) pair instability supernovae (Fig. 2). 
Note that at low metallicities, where the IMF is 
skewed toward massive stars, the fraction of mas- 
sive stars that form black holes is nearly twice as 
large as that predicted by a Salpeter IMF. Most 
of these black holes are formed through direct col- 
lapse. 

If the mass limit at which weak supernovae oc- 
cur decreases from 25 Mq down to 20 M©, the frac- 
tion of neutron stars and typical Type Hp super- 
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novae at low metallicities drops below 70 %. The 
fraction of stars that form weak Hp supernovae 
and black holes increases to compensate this de- 
crease. Table 1 summarizes the population frac- 
tions for different assumptions for the IMF and 
for the lower limit of the stellar core mass (i.e., 
lower limit of the initial mass for hydrogen-covered 
stars) resulting in fallback black hole formation. 
As mentioned above, here we assume that this cor- 
responds to the maximum stellar/core mass form- 
ing strong SNe. 

In Panel A of Fig. 7, we show the distribution of 
Type II supernovae. Most 90 %) single massive 
stars produce Type II SNe {solid line). Most of 
these produce normal Type Hp SNc {dashed line). 
Roughly 10% of all massive stars produce weak 
Type Hp SNe {dot-dashed line). As the metallic- 
ity approaches solar, some fraction of massive stars 
will produce Type IIL SNe. In Panel B, we plot 
the Type Ib/c SNe distribution. Single stars will 
not produce Type Ib/c SNe until the metallicity 
gets large enough to drive strong winds. At first, 
most Type Ib/c SNe will be produced by "weak" 
explosions that form black holes by fallback {dot- 
dashed line), but as the metallicity rises, an in- 
creasing fraction of "strong" Ib/c SNc is produced 
{long dashed line). Pair- instability SNc only oc- 
cur at low metallicities and. for our choice of IMF, 
both pulsational and non-pulsational pair instabil- 
ity supernovae each constitute only about 1 % of 
all massive stars. When using the IMF by Naka- 
mura & Umemura (2001) the pair SNe rates in- 
crease by a factor ~ 3 {thin lines). Note that in 
Panel B of Fig. 7 the pair SNe rate is scaled by a 
factor 10. 

Fig. 8 shows the distribution of GRBs and Jet- 
SNe (explosions arising from coUapsars). Since 
in the frame of the present paper we cannot well 
distinguish between GRBs and JetSNe and, lack- 
ing a better understanding of rotation, these rates 
are upper limits only. The solid line in Fig. 8 re- 
flects the total fraction of massive, single stars that 
could produce GRBs or JetSNe. The dotted line 
denotes the fraction of massive stars that could 
produce GRBs. To produce GRBs, the massive 
star must lose all of its hydrogen envelope, but 
still collapse to form a black hole. Hence, there 
is a narrow window of metallicities which allow 
GRB production in single stars. Because pulsa- 
tional instabilities arc able to eject the hydrogen 



envelope of stars even at zero metallicities, some 
GRBs could be formed at low metallicities. As in 
Fig. 6, the thin lines denote the differences caused 
by using the Nakamura & Umemura (2001) IMF 
at low metalicities. 

To determine a distribution of evolutionary out- 
comes versus redshift, we not only need to know 
the metallicity dependence of stellar winds, but 
we also need to know the metal distribution and 
spread as a function of redshift. This cosmic age- 
metallicity relation is likely to have large spreads 
and a weak trend (Pei & Fall 1995), as is also the 
case for this relation within the Milky Way (Mat- 
teucci 2001; Pagel 1997). These dependencies are 
difficult to determine because on a more global 
galactic or cosmological scale metals may be re- 
distributed so that, e.g., most of the metals even 
for low metallicity stars could be produced in stars 
of metallicity. However, to give a flavor of possi- 
ble redshift effects, we assume that the metallicity 
axis in Figs. 1-4 is indeed logarithmic and use the 
metallicity redshift distribution assumed by Lloyd- 
Ronning et al. (2002): Pei et al. (1999) distribu- 
tion versus redshift with a Gaussian spread using 
a 1 — (T deviation set to 0.5 in the logarithm of 
the metallicity. With these assumptions we can 
determine the distribution of neutron stars (thick 
solid line), black holes (thick dotted line). Type II 
SNe (thin solid hue). Type Ib/c SNe (thin dotted 
line), pair supernovae (thin dashed line) as a func- 
tion of redshift (or look-back time; Fig. 9). This 
suggests a trend in the populations of massive star 
outcomes versus redshift. 

6. Uncertainties and possible consequences 
6.1. Uncertainties in mass loss 

Our mass loss rates explicitly include only ra- 
diatively driven mass loss, though the exact na- 
ture of the Wolf Rayet star mass loss is unknown. 
We do not include pulsational ejection and similar 
eruptions or by excretion disks in rapidly rotating 
stars ("ri-limit"; Langer 1997). The magnitude 
of these mass loss mechanisms depends upon the 
composition of the star. For hot stars both the 
absolute value and the metallicity-dependencc of 
wind-driven mass loss are reasonably well under- 
stood and theoretically modeled (Kudritzki 2000, 
2002). For most of the other mass loss mecha- 
nisms and temperature and mass regimes, we have 
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insufficient observational data or theoretical mass 
loss models to make precise predictions of a mas- 
sive star's destiny. This is one reason we do not 
give precise values for mctallicity along the axes 
in Figs. 1-4. 

Though there is general consensus that reduc- 
ing the initial mctallicity of a massive star will 
increase its mass when it dies, the scaling of mass 
loss with Z during different stages of the evolu- 
tion is controversial. We have made the simplest 
possible assumption, that mass loss rates scale 
everywhere as the square root of initial metal- 
licity, essentially as the square root of the iron 
abundance. This is almost certainly naive. Vink 
et al. (2001) argue for a scaling Z°-^^ for stars with 
Toff > 25, 000 K and Z^-^^ for B-supergiants with 
Teff < 25, 000 K. Nugis & Lamcrs (2000) argued 
for a Z"-^ scaling in WN and WC stars, but for WC 
stars at least they had in mind the abundance of 
carbon in the atmosphere of the star, not the ini- 
tial mctallicity. On theoretical grounds, Kudritzki 
(2002) discusses a universal scaling for mass loss 
in hot stars that goes at Z^-^ but which has a 
threshold below which the mass loss declines more 
sharply. 

For red supergiants, even the mass loss at so- 
lar metallicity is not well determined. At higher 
stellar masses the mass loss from luminous blue 
variables and WR stars also constitutes a major 
source of uncertainty as do pulsationally-induccd 
and rotationally- induced outflows (sec above). 

6.2. Uncertainty in the explosion mecha- 
nism 

The mechanism whereby the collapse of the 
iron core in a massive star results in a strong 
explosion has been debated for decades. The 
current paradigm is based on a neutrino pow- 
ered "hot bubble" formed just outside the young 
proto-neutron star, but even the validity of this 
paradigm is debated along with its specific pre- 
dictions (Herant et al. 1994; Burrows et al. 1995; 
Janka & Mueller 1996; Mezzacappa et al. 1998). 
The role of rotation and magnetic fields is also con- 
tentious (Leblanc & Wilson 1970; Fryer & Heger 
2000; Ardeljan et al. 2001; Wheeler et al. 2002; 
§ 6.3). 

Our intuition here has been guided by para- 
metric surveys in which the explosion is simu- 



lated using a piston. The numerous uncertain- 
ties are thus mapped into choices of the piston's 
location and motion. These parameters arc con- 
strained by the requirement that the explosion 
not eject too much neutron-rich material (hence 
a minimum mass interior to the piston) and that 
the kinetic energy of the explosion measured at 
infinity be 10^^ erg. Though a single event, SN 
1987A occurred for a representative helium core 
mass (6Mq) and had a measured kinetic energy 
at infinity of ~ 1 - 1.5x10^1 erg (Woosley 1988; 
Arnett et al. 1989; Bethe & Pizzochero 1990). The 
requirement that supernovae typically make ^0.1 
M0 of ^^Ni also means that the piston cannot be 
situated too far out or produce too weak an explo- 
sion. There are also more subtle conditions - that 
the mass cut frequently occur in a location where 
past (successful) calculations of the explosion have 
found it, that the distribution of remnant masses 
resemble what is observed for neutron stars, that 
the integrated ensemble of abundances resemble 
Population I in our galaxy, and so on. 

Fig. 5 shows the remnant masses for a survey 
of explosions in solar metallicity stars that ne- 
glects mass loss. The progenitor stars described 
in Woosley et al. (2002) were exploded using a 
piston located at the edge of the "iron core" . The 
iron core was defined by the location of an abrupt 
jump in the neutron excess (electron mole number 
= = 0.49). A constant kinetic energy at infinity 
(1.2x 10^^ erg) was assumed (see also Woosley & 
Weaver 1995). In fact, the explosion energy will 
probably vary with mass. Fryer (1999) calculates 
that the explosion energy will actually weaken as 
the mass of the helium core increases. Thus fall 
back could have an even earlier onset and more 
dramatic effects than Fig. 5 would suggest. 

The apparent non-monotonic behavior in Fig. 5 
is largely a consequence of the choice of where the 
piston was sited. The neutronized iron core may 
have a variable mass that depends on details of 
oxygen and silicon shell burning (Woosley et al. 
2002). The density gradient around that core can 
also be highly variable. Thus enforcing a constant 
kinetic energy at infinity does not always lead to 
a predictable variation of remnant mass with ini- 
tial mass. More recent unpublished calculations 
by Heger et al. (2003a), also of zero metallicity 
stars, place the piston at an entropy jump (di- 
mensionless entropy S/N^kB ~ 4) rather than a 
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% jump. This choice, which is more consistent 
with explosion models, assumes that an explosion 
develops when the accretion rate declines rapidly. 
The rapid decline is associated with the density 
(and entropy) discontinuity near the base of the 
oxygen burning shell. Such a prescription gives 
more nearly monotonic results and, in particular 
the bump around 17 M0 in Fig. 5 is absent. 

Nevertheless Fig. 5 does suggest that the lines 
separating black hole formation by fall back from 
neutron stars in Figs. 1-4 should be interpreted 
only as indicating trends. They may not be as 
smooth or as monotonic as indicated. 

6.3. Uncertainty in the effects of rotation 

Rotation can enhance the mass loss in stars and 
a spread in initial rotation can smear out the tran- 
sitions between the different mass and metallic- 
ity regimes. We have not considered cases where 
rotationally enhanced mass loss might be impor- 
tant. In such cases the limiting mass for loss of the 
hydrogen envelope could be lowered and, at the 
same time, the mass of the helium core increased 
(Heger et al. 2000; Meynet & Maeder 2000). The 
higher mass loss would tend to lower the metallic- 
ity for divisions between Type II and Type Ib/c 
supernovae as well is the divisions between strong, 
weak and no supernova explosions. The higher he- 
lium core masses with increase the metallicity di- 
visions between strong, weak, and no supernova 
explosions. The total change will depend on the 
competition of the larger helium core masses and 
enhanced mass loss rate. 

If the core is rotating rapidly at collapse, ro- 
tation may also influence the explosion mecha- 
nism and especially the possibility of making a 
GRB. Also pair-creation supernovae could be sig- 
nificantly affected by rotation, in particular the 
lower mass limit for direct black hole formation 
(Glatzel et al. 1985; Stringfellow & Woosley 1988). 
Early calculations that followed angular momen- 
tum in massive stars (e.g., Kippenhahn & Thomas 
1970; Kippenhahn et al. 1970; Endal & Sofia 1976, 
1978; Tassoul 2000) all found sufficient angular 
momentum retained in the core to reach criti- 
cal rotation ("break-up velocity") before the fi- 
nal central burning phases. More recent calcu- 
lations by Heger (1998); Heger et al. (2000) find 
presupernova core rotation rates in massive stars 
that would lead to sub-millisecond neutrons stars 



just around break-up if angular momentum were 
conserved perfectly during the collapse. Calcula- 
tions by Meynet & Maeder (1997); Maeder & Zahn 
(1998); Maeder & Meynet, priv. com. (2000) indi- 
cate core rotation rates after central helium burn- 
ing similar to those found by Heger et al. (2000). 

Recently Spruit (2002) has discussed a "dy- 
namo" mechanism based on the interchange in- 
stability that allows the estimation of magnetic 
torques to be included in models for stellar evo- 
lution. Preliminary calculations by Heger et al. 
(2002b, 2003b) that include these torques find a 
presupernova angular momentum equivalent to 5 
-10 milliseconds - still somewhat faster than ob- 
served young pulsars, but too slow for coUapsars. 
If the estimates of magnetic torques by Spruit 
(2002) are valid then single stars are unlikely to 
produce coUapsars and rotation is probably not 
a factor in the explosion of common supernovae. 
Nevertheless, in Fig. 3 we indicate the regimes 
where the structure of the star, excluding the ques- 
tion of sufRcient rotation, is favorable for coUap- 
sars and GRBs. 

7. Conclusions and Observational Tests 

We have described, qualitatively, the likely fates 
of single massive stars as a function of metallicity. 
Our results suggest various trends in the obser- 
vations of these objects which may be subject to 
observational tests. 

• Normal Type Ib/c SNe are not produced by 
single stars until the metallicity is well above 
solar. Otherwise the helium core mass at 
death is too large. This implies that most 
Type Ib/c SNe are produced in binary sys- 
tems where the binary companion aids in re- 
moving the hydrogen envelope of the collaps- 
ing star. 

• Although less extreme than Type Ib/c SNe, 
single stars also do not produce Type IIL 
SNe at low metallicities. Similar to Type 
Ib/c SNe, Type IIL SNe from single stars are 
probably "weak" SNe until the metallicity 
exceeds solar, also implying that Type IIL 
SNe are produced in binaries. 

• If GRBs are produced by single star collapse 
(perhaps unlikely given the constraints on 
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angular momentum), single stars only make 
up a small subset of GRB progenitors at 
higher metallicitics. It is more likely that 
binary systems form GRBs. Such systems 
will occur more frequently at low metallici- 
ties (Fryer et al. 1999). 

• Jet-driven supernovae from single stars are 
likely to be much more common than GRBs 
from single stars. 

It is difhcult to make direct comparisons to ob- 
servations without including binary stars in our 
analysis, but there are a number of constraints 
that should be considered. First, an increasing 
number of JetSNe and weak supernovae explosions 
are being discovered (Nakamura et al. 2001; SoUer- 
man et al. 1998; Turatto et al. 1998). Although 
there is an observational bias against the discov- 
ery of weak supernovae and they are much dimmer 
than JetSNe, they may still dominate the sample 
of stars more massive than 25 M©. Clearly, good 
statistics (and correct analysis of the systematics) 
are necessary to determine the relative ratio of jet- 
driven and weak SNe. With such statistics, we 
may be able to place constraints on the rotation 
of massive stellar cores. 

If the IMF becomes more top-heavy at low 
metallicity (< 10"'' Z©; Bromm et al. 2001; Schnei- 
der et al. 2002) the number of core collapse super- 
novae (mostly Type Hp) and GRBs (if occurring 
in single stars) should significantly increase at 
high redshift. If the current estimates of a char- 
acteristic mass of ~ 100 M© for primordial stars 
(Bromm et al. 1999; Abel et al. 2002; Nakamura 
& Umemura 2001) is correct we should expect a 
large fraction of pair SNe and very massive black 
holes (or Type III collapsars) at zero metallicity, 
as well as an increase of massive black holes from 
stars in the 60-140 Mq region. 
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Table 1: Remnant and supernova population yields for different metallicitics, IMFs, and mass limits. 



Zero Mctallicity Solar Metallicity 

Object _higliMi,™H— 1owA4™h high A4™h low il4™H 

IMFsal IMFnu IMFsai IMFnu IMFgal IMFsal 



Remnants 


NS 


75 


56 


66 


50 


87 


75 


BH 


23 


36 


32 


43 


13 


25 


MBH 


0.9 


3.0 


0.9 


3.0 








Supernovae 


lip Strong 


75 


56 


66 


50 


77 


70 


Hp Weak 


12 


8.9 


21 


16 





6.9 


IIL 














6.4 


6.4 


Ib/c Strong 














9.2 


5.1 


Ib/c Weak 














7.6 


12 


Other Outbursts 


Puis. Pair 


1.4 


4.7 


1.4 


4.7 








Pair SNe 


1.4 


4.6 


1.4 


4.6 








Jet SNe 


24 


39 


33 


46 


13 


25 


GRBs 


1.4 


3.4 


1.4 


4.7 


7.8 


12 



Note: For solar metallicity we use the IMF by Salpeter (1955; IMFgai), for zero metallicity we additionally 
supply the results for the IMF by Nakamura & Umemura (2001; IMFnu)- We give the results two different 
lower mass limits for fallback black hole formation (A/pgjj): high corresponds to 25 Mq and low to 20 Mq 
(Fryer 1999). 
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Fig. 1. — Remnants of massive single stars as a function of initial mctallicity {y-axis; qualitatively) and 
initial mass (x-axis) . The thick green line separates the regimes where the stars keep their hydrogen envelope 
(left and lower right) from those where the hydrogen envelope is lost (upper right and small strip at the 
bottom between 100 and 140 M0). The dashed blue line indicates the border of the regime of direct black 
hole formation (black). This domain is interrupted by a strip of pair-instability supernovae that leave no 
remnant (white). Outside the direct black hole regime, at lower mass and higher mctallicity, follows the 
regime of BH formation by fallback (red cross hatching and bordered by a black dash-dotted line). Outside 
of this, green cross hatching indicates the formation of neutron stars. The lowest-mass neutron stars may 
be made by 0/Ne/Mg core collapse instead of iron core collapse (vertical dash-dotted lines at the left). At 
even lower mass, the cores do not collapse and only white dwarfs are made (white strip at the very left). 
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Fig. 2. — Supcrnovae types of non-rotating massive single stars as a function of initial mctallicity and initial 
mass. The lines have the same meaning as in Fig. 1. Green horizontal hatching indicates the domain where 
Type Hp supernovae occur. At the high-mass end of the regime they may be weak and observationally faint 
due to fallback of ^^Ni. These weak SN Type Hp should preferentially occur at low metallicity. At the 
upper right edge of the SN Type II regime, close to the green line of loss of the hydrogen envelope, Type 
IIL/b supernovae that have a hydrogen envelope of < 2M0 are made (purple cross hatching). In the upper 
right quarter of the figure, above both the lines of hydrogen envelope loss and direct black hole formation. 
Type Ib/c supernovae occur; in the lower part of their regime (middle of the right half of the figure) they 
may be weak and observationally faint due to fallback of ^^Ni, similar to the weak Type Hp SNe. In the 
direct black hole regime no "normal" (non-jet powered) supernovae occur since no SN shock is launched. 
An exception are pulsational pair-instability supernovae (lower right corner; brown diagonal hatching) that 
launch their ejection before the core collapses. Below and to the right of this we find the (non-pulsational) 
pair-instability supernovae (red cross hatching), making no remnant, and finally another domain where black 
hole are formed promptly at the lowest metallicities and highest masses (while) where nor SNe are made. 
White dwarfs also do not make supernovae (white strip at the very left). 
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9 10 25 34 40 60 100 140 260 

initial mass (soiar masses) 

Fig. 3. — CoUapsar types resulting from single massive stars as a function of initial metallicity and initial 
mass. Lines have the same meaning as in Fig. 1. Our main distinction is between coUapsars that form from 
fallback (Type II; red) and directly (Type I; pink). We subdivide these into those that have a hydrogen 
envelope {cross hatching), only able to form jct-powcrcd supcrnovac (JetSNe) and hydrogen-free coUapsars 
{diagonal cross hatching), possibly making either JetSNe or GRBs (see also Fig. 4). The first subclass is 
located below the thick green hne of loss of the hydrogen envelope and the second is above it. The hght 
brown diagonal hatching at high mass and low metallicity indicates the regime of very massive black holes 
formed directly (Type III coUapsars) that collapse on the pair-instability and photo-disintegration. Since 
the coUapsars scenario require the formation of a BH, at low mass (left in the figure) or high metallicity (top 
of the figure) and in the strip of pair-instability supernovae (lower right) no coUapsars occur {white). 
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Fig. 4. — Jet-driven supernovae types as a function of initial mctallicity and initial mass. Lines have the 
same meaning as in Fig. 1. The regimes in which hydrogen-rich JetSNe are possible (below the thick green 
line indicating loss of the hydrogen envelope) is indicated by cyan hatching, and that of hydrogen- free JetSNe 
by Ught brown hatching (above the thick green Une) . In the latter regime also GRBs may be possible, while 
in the first regime a hydrogen envelope is present and the travel time of a relativistic jet though it is much 
bigger than typical observed GRB durations. In the region of very massive black hole formation {magenta 
cross hatching; lower right corner) long JetSNe and long X-ray outbursts may occur since the bigger mass- 
scale of these objects also translates into a longer time-scale. If these objects are at cosmological distances, 
additionally the apparent time-scale and wavelength are both is stretched. 
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Fig. 5. — Remnant masses of metal- free stars as a function of initial mass for stars from Woosley et al. 
(2002, WHW02, solid line) assuming a constant kinetic energy of the ejecta of 1.2x10^^ erg. The explosions 
were simulated by a piston at the edge of the deleptonized core similar to Woosley & Weaver (1995, WW95, 
crosses) . 
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Fig. 6. — Fraction of massive stars that form neutron stars {solid line) and black holes {dotted line) as a 
function of metallicity for a Salpeter initial mass function (thick lines; Salpeter 1955). The dashed lines 
denote just those black holes formed through fallback and the dot-dashed lines denote black holes formed 
from very massive (> 300 M©) stars. The thin lines arise from assuming the IMF at low metallicities is given 
by Nakamura & Umemura (2001) at low metallicity (see § 5). Note that at low metallicities, pair-instability 
supernovae leave no compact remnant whatsoever, so that in this regime the total of all fractions is less than 
one. 
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Fig. 7. — Fraction of massive stars that form supernovae for a Salpeter initial mass function (thick lines; 
Salpeter 1955). At low redshifts we use an alternate initial mass function {thin hnes) from Nakamura & 
Umemura (2001). Most single stars become Type II supernovae {soUd hne, Panel A), and most of these are 
strong Hp SNe {dotted hne). Roughly 15 % of Type II SNe are weak Type Hp supernovae {dot-dashed hne). 
As the metallicity approaches solar, the fraction of weak supernovae decreases and a small fraction of Type 
IIL SNe are produced {dashed hne). Type Ib/c supernovae are not produced until the metallicity approaches 
solar {sohd hne, Panel B), and most of these SNe will be weak {dot-dashed hne). Not until the metallicity 
exceeds solar are strong Ib/c SNe produced {long dashed hne). Pair- instability supernovae {dashed line, 
Panel B) and pulsational pair-instability supernovae {dotted line, Panel B) are rare and only produced at 
low metallicitics. Their fraction depends strongly on the unknown IMF at these low metallicitics. Note that 
in the figure wc multiply the pair instability SNe fraction by a factor 10. 
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Fig. 8. — Upper limits on the fraction of massive stars that form jet-driven supernovae and gamma-ray bursts 
for a Salpcter initial mass function {thick lines; Salpeter 1955). At low rcdshifts we use an alternate initial 
mass function {thin lines) from Nakamura & Umemura (2001). These upper limits are determined assuming 
all massive stars have the necessary rotation rates to produce coUapsars. Single stars produce GRBs mostly 
in a narrow range of metallicities, but can produce Jet SNe at all metallicities until the metallicity is so high 
that mass loss prohibits the formation of black holes. 
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Fig. 9. — The distribution of neutron stars {thick solid line), black holes {thick dotted line), Type II SNe 
(thin solid line), Type Ib/c SNe {thin dotted line), pair supernovae {thin dashed line) as a function of redshift. 
We have assumed that the metallicity axis in Figs. 1 and 2 is indeed logarithmic with the maximum mass 
for which pair creation supernovae occur (Fig. 2) corresponding to a metallicity of lO^** solar. We have used 
the metallicity redshift distribution assumed by Lloyd-Ronning et al. (2002): Pei et al. (1999) distribution 
versus redshift with a Gaussian spread using a 1 — ct deviation set to 0.5 in the logarithm of the metallicity. 
This gives an idea of the trends in the populations of massive single star outcomes as a function of redshift. 
Given the various assumptions that have to be made for this type of analysis, these the resulting absolute 
numbers should be interpreted with a great caution. 
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